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Abstract 

Many-multiplet (MM) method applied to three inhomogeneous 
samples of Keck/HIRES quasar absorption spectra gives a shift in the 
value of the fine-structure constant of Aa/a = (—5.4 ± 1.2) x 1CP 6 in 
the redshift range 0.2 < z < 3.7 The 1 a error claimed in [2] is, 
however, much too small and cannot be maintained by current obser- 
vations of quasars. We present a modified MM method to set an upper 
limit on Aa/a from a homogeneous sample of Fen lines identified in 
the up-to-date best quality VLT/UVES spectrum of HE 0515-4414. 
Our result is Aa/a = (1.1± 1.1) x 10 -5 at z = 1.149. Theoretical mod- 
els of the fundamental physical interactions predict that the proton- 
to-electron mass ratio (// = m p /m e ) may relate to the shift in Aa/a 
as Afj,/ fi = TZAa/a. We use VLT/UVES high-resolution observations 
of molecular hydrogen H2 ultraviolet absorption lines at z = 3.025 to- 
ward Q 0347-3818 to bound the value of 1Z. The obtained constraints 
on A In \i = (2.1 ± 3.6) x 10~ 5 and on A In a = (1.1 ± 1.1) x 10~ 5 rule 
out vary large values of \1Z\ > 6. Future observations with a new High 
Accuracy Radial velocity Planet Searcher (HARPS) spectrograph may 
provide a crucial test for the Aa/a measurements at a level of 10~ 6 . 

1 Introduction 

The purpose of this paper is to outline some problems related to the mea- 
surements of the dimensionless fundamental constants from quasar (QSO) 
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absorption-line spectra: the fine-structure constant a = e 2 /Tic and the proton- 
to-electron mass ratio \i = m p /m e , which characterize the strength of elec- 
tromagnetic and strong interactions, respectively. 

The current interest in this field stems from the prospect of using the 
high-resolution spectra of extragalactic objects to measure with high accu- 
racy the transition energies of a certain ionic species and molecules at differ- 
ent cosmological redshifts z 1 . If physical constants varied in the past, then 
the comparison of the transition energies measured from QSO spectra with 
their laboratory values provides direct astrophysical constraints on the ra- 
tios a z /a and fi z / fi over a cosmological time-scale (At ~ 10 10 yr). We define 
Aa/a = (a z — a) /a and Afij fi = (fi z — fi)/fi, where index z marks the 
values of the fine-structure constant and the proton-to-electron mass ratio 
at redshift z, whereas their laboratory values at z = are denoted by a and 
/i, respectively 2 . At present, these constants are measured with the relative 
accuracies 5 a = 4 x 1CT 9 and 5^ = 2 x 10~ 9 and their values are equal to 
a = 1/137.03599958 and fi = 1836.1526670 0. 

The variability of the fundamental constants was firstly suggested by 
Milne |3] and Dirac [5] who assumed that the Newtonian gravitational con- 
stant G may change in cosmic time. Different functional forms for the time 
variations of G and other constants were later considered in a series of pub- 
lications (see 0, for a review). The variability of a caused by a coupling 
between electromagnetism and gravity is considered within the framework of 
different modifications of the Kaluza-Klein theory (e.g., 0-JT2])- The num- 
ber of theoretical publications on this topic is considerably increased in the 
last two years. 

For astronomical observations, it is more convenient to deal with dimen- 
sionless constants since they are independent on the choice of physical units. 
That is why most efforts to reveal hypothetical variability of constants are 
connected with measurements of a (see and references cited therein). 

First observational study (1967) of the time dependence of a from QSO 
absorption spectra was based on the analysis of alkali doublets (AD method) 
and yielded Aa/a = —0.02 ± 0.05 at z = 1.95 [Hj. Since then, there have 
been made many other AD studies of Aa/a with increasing accuracy in 

1 z = (A — Ao)/Ao, where A is the observed wavelength of a given ionic transition and Ao 
is the corresponding rest frame (laboratory) value. Modern observational data are ranging 
from z = (the Milky Way) up to z ~ 6 (distant QSOs). 

2 In the literature fi is sometimes referred to its inverse value, /i* = m c /m p . In this 
case A[i*/[i* = —A^i/^i. 



2 



accord with increasing quality of observations: Aa/a = (0 ± 1.5) x 10 3 at 
z = 3.2 Aa/a = (-4.5 ± 4.3 ± 1.4) x 10" 5 at z = 2 - 4 [Sj, Aa/a = 
(—0.5 ± 1.3) x 10 -5 at z = 2 — 3 17J. However, the use of absorption lines 
from different multiplets in different ions, the so-called many-multiplet (MM) 
method, may in principle increase the accuracy of these Aa/a estimations 



Applying the MM method to three samples of Keck/HIRES QSO data, 
Murphy et al. [2| obtained a statistically significant negative value of Aa/a = 
(-5.4 ± 1.2) x 10~ 6 in the redshift interval 0.2 < z < 3.7. In spite of all ad- 
vantages of the MM method, the result is nevertheless surprising since it im- 
plicitly assumes extremely stable conditions during the observing nights and 
unrealistically homogeneous physical parameters in the absorption clouds. 
We consider computational caveats of the AD and MM methods in Section 2. 

It should be emphasized that recent improvements in meteoritic abun- 
dance determinations lead to the bound Aa/a < 3 x 10 -7 at z — 0.45 [22] 
which contradicts the MM measurements [2] at the same redshift. 

If a is supposed to be time dependent, the other gauge coupling constants 
should also depend on time as predicted by theoretical models of the funda- 
mental physical interactions (e.g., JE]-[2II])- However, the relation between 
the variations of the low energy gauge couplings is highly model dependent. 
For instance, within the framework of a unified theory [e.g., supersymmetric 
SU(5)] the time variation of the QCD scale parameter A leads to the change 
of the proton mass m p which is proportional to the shift Aa/a [T%1 FEU l2Tij: 



where the factor 1Z is known with large theoretical uncertainties |2*3^ 12*4^ 12*5] . 

Theoretically predicted numerical value of 1Z can be tested by measuring 
the electronic-vibrational-rotational lines of the H2 molecule observed in QSO 
spectra at high- ,2 as described in j2H]- We discuss this test in Section 3. Our 
conclusions and remarks for future prospects are summarized in Section 4. 
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2 Methods to constrain Aa/a from QSO ab- 
sorption spectra 



2.1 The alkali-doublet (AD) method 



The value of the fine-structure constant a z at redshift z can be directly 
estimated from the relative wavelength separation between the fine-splitting 
lines of an alkali doublet (AD method) [27| : 



where AX Z and (X) z are, respectively, the fine-structure separation and the 
weighted mean wavelength for a given doublet from a QSO absorption system 
at redshift z, while AA and (A) denote the same quantities from laboratory 
data. 



where 5\ z = <T\ z /A\ z , 5\ = a\/A\, and a\ z , a\ are, respectively, the errors 
of the wavelengths in the absorption system and laboratory. 

Equation (J2J) shows that the AD method uses only the wavelength ratios 
and thus it does not require the measurements of absolute wavelengths. This 
fact excludes large systematic errors caused by the wavelength calibration in 
a wide spectral range. Therefore the accuracy of the AD measurements is 
mainly restricted by statistical fluctuations in the recorded counts at different 
wavelengths within the line profile. 

The measurement of the line profile parameters and their errors are thor- 
oughly described in [21]. The midpoint A c of an interval over M pixels 
covering the line profile can be calculated from the equations 





The uncertainty a& a / a in the individual measurement can be estimated from [2H] : 





Ac = ™r£Mi-va) 



(4) 



and 



M 



W x = a 



(5) 
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where W\ is the equivalent width, a is the wavelength interval between pixels, 
Ii is the observed count rate at pixel i, is the continuum count rate at 
pixel i, and Aj is the wavelength at pixel i. 

The errors in A c and W\ caused by counting statistic and the uncertainty 
of the overall height of the continuum level are given by 



0% 



a 



%) 2 [A - 2X C B + MA C 2 ) + (^) 2 (B - MA C ) 2 



.1/2 



(6) 



and 



(aj/Cf M + (a c /C) 2 (M - W x /a) 2 } 1/2 , (7) 



where A = X^Af, B — (&i/C) is the inverse mean signal-to-noise 

(S/N) -1 ratio over the line measurement interval, and {ac/C) is the mean 
accuracy of the continuum level. 

If the rms error of the wavelength scale calibration is cr sys , then the total 
error is given by the sum 

<t = < + <s- (8) 

The error of the wavelength scale calibration is related to the radial velocity 
accuracy as 

a v = ca sys /X , (9) 

where c is the speed of light. 

For example, a sys ~ 0.002 A for VLT/UVES data |2D] corresponds at 
5000 A to cr sys = 0.12 km s -1 . If a Xc = 0.002 A, then for the Mgll 
AA2796,2803 A doublet © gives a Aa/a = 1.3 x 10~ 4 at z = 2, which is 
a typical accuracy for the individual AD estimations from isolated doublets 
with W\ ~ 100 mA. Equation © shows that a\ c decreases linearly with W\, 
meaning that, not surprisingly, the error for weak absorption lines (W\ 1 
A) may be very large. For example, the rest-frame equivalent widths for Fell 
A1608 A and A1611 A at z = 3.39 toward Q0000-2620 measured from the 
Keck/HIRES data are W W08 ~ 0.178 A and W 1611 ~ 0.003 A [3T|, i.e. the 
measured errors a\ c for these lines should differ by ~ 50 times. It is clear that 
to achieve the highest accuracy in the Act /a estimations by means of the the 
AD method only well separated and unsaturated doublets with W\ ~ 100 
mA should be selected for the analysis. 
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2.2 The many-mult iplet (MM) method 

The many-multiplet (MM) method is a generalization of the AD method 
(see [2] and references therein). The method uses many transitions from 
different multiplets and different ions observed in a QSO absorption system. 
The line profiles selected for the MM analysis are represented by a sum of 
the same number of components (Voigt profiles) each of which is defined 
by three individual parameters (the column density, the Doppler width or 
6-parameter, and the redshift Zi) and a common value of Aa/a. The rest 
wavenumbers in the MM method can be altered according to [2] 

uj z =u + qx z , (10) 

where 

x z = (a z /a) 2 - 1 , (11) 

and lu z , Uq are the rest-frame wavenumbers of a transition in a QSO system 
with redshift z and measured in the laboratory, respectively. The sensitivity 
coefficients q (determining the sensitivity of u to the variation of a) are 
listed in Table 2 in 0. 

All parameters are estimated through the minimization of x 2 simultane- 
ously for the whole set of lines. A significant improvement in the accuracy 
of the Aa/a estimations by the MM method compared to the AD method 
is due to a wide range of the q values (see Section 2.3 below). Nevertheless, 
the MM method has some immanent shortcomings which could affect the 
claimed accuracy cr/Aa/a) — 1-2 x 10~ 6 . 

First of all, the weakest point of the MM method is the multiple velocity 
component Voigt profile deconvolution. From mathematical point of view, 
any deconvolution itself is a typical ill-posed problem When several 

smoothing operators (like convolution with the spectrograph point-spread 
function, summing, and integrating) are involved in the minimization proce- 
dure, the ill-posing may cause quite ambiguous results. In case of complex 
line profiles when many components are required to describe the whole pro- 
file, the fitting parameters cannot be estimated with high accuracy because 
of strong inter-parameter correlations. For example, the value of Aa / a can 
be artificially constrained by the increase of the component number or by 
variations of 6-values. Thus, to enhance the reliability of the results only 
simple line profiles (requiring minimum number of components) should be 
taken into MM analysis. 
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From physical point of view, the Voigt profile deconvolution assumes that 
the complex profile is caused by separate clouds lined along the line of sight 
with different radial velocities, each having its own constant gas density, 
kinetic temperature and, hence, its own ionization state. However, obser- 
vations of QSOs have shown that with increasing spectral resolution more 
and more subcomponents appear in the line profile. This complexity indi- 
cates that, in general, metal lines arise in continuous turbulent media with 
fluctuating velocity and density fields. 

We note in passing that the highest spectral resolution achieved in obser- 
vations of QSOs at modern giant telescopes is ~ 5 km s _1 (FWHM), whereas 
the expected values for minimum (thermal) widths of metal lines are less or 
about 1 km s" 1 for typical kinetic temperatures 100 K < T^ in < 10 4 K. This 
means that we are still not able to measure directly the true intensity in the 
line profile but observe only the apparent intensity which is a convolution of 
the true spectrum and the spectrograph point-spread function. 

Thus, to measure the transition energies, the distribution function of the 
radial velocities within the intervening cloud must be known for the AD 
method, whereas the MM method requires in addition the knowledge of the 
density distribution, for it deals with ions of different elements which may 
not arise co-spatially (SHI- Since the velocity and density distributions are 
determined from the same line profiles which are used in the measurements of 
Act /a, it is not possible to discriminate among kinematic effects and possible 
changes in a on the basis of only one absorption system. But if we have many 
absorption systems, statistical approach may probably solve this problem if 
shifts in the measured Aa/a values are of random nature and, thus, can be 
averaged out in a large ensemble of absorption systems where any deviations 
in velocity distributions of different ions are equally probable. 

Shifts in Aa / a may also arise if the main assumption of the MM method 
that different ionic species trace exactly each other in their space distribution 
is not adequate. Dissimilarities in the line profiles of different ions may be 
caused by density fluctuations leading to inhomogeneous photoionization in 
the absorption cloud. 

The QSO absorption clouds is a very inhomogeneous population. Their 
linear sizes, L, can be estimated from the observed neutral hydrogen column 
densities, iV(Hi) ~ 10 20 — 10 21 cm -2 , and the estimated gas number densities, 
n H ~ 10~ 2 — 1 cm -3 , giving L = A^(Hl)/n H ~ 10s pc - 10s kpc. The 
intervening clouds may also contain small dense filaments with much higher 
densities (nn ~ 10 2 cm -3 ) and the line-of-sight sizes less than 1 pc [3%j . 
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It is usually assumed that the metal absorbers in QSO spectra are in 
photoionization equilibrium with the ambient metagalactic UV background. 
To describe the thermal and the ionization state of the gas in the diffuse 
clouds it is convenient to introduce the dimensionless 'ionization parameter' 
U = Tiph/^H, equal to the ratio of the number of photons in 1 cm 3 with 
energies above one Rydberg to the gas number density 

For illustrative purpose, we will assume that the UV background is not 
time or space dependent and that its shape and intensity is given by the 
metagalactic spectrum at z = 3 presented in [33] (n p h ~ 2 x 10~ 5 cm -3 ). 
Then fractional ionizations T 0i j = n a ^/n a (with n a being the local number 
density of element 'a' and n a ^ the density of ions in the ith ionization state) 
of different elements can be calculated as functions of U using the photoion- 
ization code CLOUDY gE]. 

Figure 1 shows these results for a metallicity Z = 0.1Z & (typical for the 
QSO metal absorbers) and for the elements involved in the MM analysis 
in Similar profiles for different ions can only be observed when the 
ratios T^i/T^j are constant, i.e., the curves in Fig. 1 are parallel. This 
condition is realized for Mgl, Mgll, Alii, Sill, Fell, Nill, and Znll in the 
range log £7 < —4.2, i.e., for the most dense volumes of the cloud. Such dense 
clumps (filaments) have, however, a negligible filling factor since their sizes 
are small, and, thus, the most of the diffuse cloud is filled up by a rarefied, low 
density gas characterized by higher values of log [7. In the range log £7 > —4, 
the diversity of the Tx values is considerable and one may expect to observe 
nonsimilar profiles of the above mentioned ions. Under these circumstances, 
in order to be more secure the MM method should work with homogeneous 
samples of ions which have undoubtly the same volume distribution like, e.g., 
Fe II lines. 

Note that Aim is obviously not a good candidate for MM calculations 
since the shape of the Taihi curve clearly differs from the others. The fact 
that no additional scatter in Aa/a for systems containing Aim was found 
in j2] may imply that the available spectral resolution is not high enough to 
obtain the true profiles of the metal lines. 

To conclude this section we comment on the accuracy cr^ a / a which can 
be in principle achieved from the up-to-date QSO observations. Using the 
method of error propagation, one finds from (jlUj) the relative error of uj z \ 

<L, = <*wJUz = 2 \Q\ VAa/a ■ (12) 

Here Q = q/uj® is the dimensionless sensitivity coefficient. 
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Figure 1: Fractional ionizations logT^ vs. the ionization parameter \ogU 
in case of the Haardt-Madau [35] metagalactic UV background at z = 3 and 
metallicity Z = 0.1Z Q . Similarity of the line profiles of a pair of ions 'a' and 
requires the fractional ionization ratio T a /Tf, to be constant for different 

U 

For the lines listed in Table 2 in the Q values are ranging from —0.027 
(Cm A2066 A) to 0.050 (Znll A2026 A), however, for the most frequently 
used transitions in the MM calculations (Table 1 in [2 ), this range is nar- 
rower: Q(Mgll A2803 A) = 0.0034, and Q(Fell A2600 A) = 0.035. These 
numbers show that the error of the sample mean cr<Aa/a) = 1-2 X 10~ 6 found 
in j2] requires the relative error (<L Z ) — 8.6 x 10~ 9 — 8.4 x 10~ 8 , meaning that 
the sample mean wavelength accuracy should be {<j\) ~ 2 x 10~ 5 — 2 x 10~ 4 
A at 2000 A (cta = A S Uz ), which is equal to the radial velocity accuracy of 
Av ~ 2.5 — 25 m s _1 (Av = c5 Wz ). Even if the individual measurement is an 
order of magnitude less accurate as compared to the sample mean, the un- 
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Figure 2: Fell lines associated with the z = 1.149 absorption system toward 
the quasar HE 0515-4414 (normalized intensities are shown by dots). The 
mean signal-to-noise (S/N) ratios per resolution element are indicated 

certainty of Av ~ 25 — 250 m s _1 seems to be hardly realized for all systems 
collected in Namely, the thermal shifts of order 300 m s _1 are expected 
if temperature changes on 0.5 K between the science and calibration spectra 
exposures |3Uj . 

2.3 The regression MM method 

The standard MM method based on the multiple velocity component Voigt 
profile fit to the complex spectra when Aa/a is treated as an additional 
free parameter [2] can be easily modified to be free from the uncertainties 
mentioned in the previous section. By its very nature, the MM method is 
similar to that developed in for the A/x//x calculations. So, by analogy, 
we can re- write (jlOj) in the linear approximation (| Act/a | <C 1) in the form 

Zi = z + K a Qi, (13) 

where the slope parameter K a is given by 

x Aa* . s 

Ka = 2 1 + z ) . 14 

a 
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Table 1: Fell lines at z = 1.149 toward the quasar HE 0515-4414 and sensi- 
tivity coefficients Q 



Ag, A 


/ 


A C) A 


W x , mA 


z 


Q 


1608.45085(8) 


0.0580 b 


3456.7176(46) 


7.8(1.8) 


1.1490975(28) 


-0.0193(48) 


2344.2130(1) 


0.1140 6 


5037.9398(23) 


27.5(1.6) 


1.1490964(10) 


0.0294(35) 


2374.4603(1) 


0.0313 6 


5102.9319(79) 


8.4(1.5) 


1.1490913(33) 


0.0389(35) 


2382.7642(1) 


0.3006 c 


5120.8015(22) 


61.0(1.9) 


1.1491012(9) 


0.0357(36) 


2586.6496(1) 


0.0684 d 


5558.9439(29) 


20.1(1.4) 


1.1490904(11) 


0.0393(39) 


2600.1725(1) 


0.2239 c 


5588.0269(21) 


59.1(2.0) 


1.1490985(8) 


0.0353(39) 



"Vacuum rest wavelengths are taken from j2]. 



REFERENCES: 6 Welty et al. [40,; c Morton [H]; d Cardelli & Savage 



Here Aa*/a = —Aa/a. 

If Aa/a is not zero, there should be a correlation between Z{ (= \ c>i / A 0i j — 
1) and Qi, and the regression parameters Zq and be estimated from 

the measured centers of the isolated metal lines. 

For instance, Fell lines listed in Table 2 in for which the Q values 
are ranging in the broad interval between —0.019 and 0.035, may be suitable 
for such calculations. Equation (|T3"j) shows that the broader the interval of 
Q- values, the higher the accuracy of the slope parameter n a (see, e.g., J3Z] 
where the errors of the linear regression are calculated) . 

It is also an advantage that Fe II lines arise mainly in low ionized regions, 
where their thermal width is the smallest among other elements (Mg, Al, 
Si) proposed for the MM method (we do not consider Cm, Nill, and Znll 
as appropriate for the MM analysis because they are rare in QSO spectra 
and, when observed, are usually very week) since the narrow lines allow to 
measure the transition energies with highest accuracy. Below we consider an 
example of such Fe II sample. 

Figure 2 shows isolated and unsaturated Fell lines from the metal sys- 
tem at z — 1.149 toward the bright quasar HE 0515-4414. The spectra of 
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Figure 3: Relation between the redshifts z\ calculated for individual Fell 
lines listed in Table 1 and their sensitivity coefficients Qi. The solid line is 
the linear regression z f = z + n a Qi- The dotted lines representing the 1 o 
deviations from the optimal slope demonstrate that Aa/a = at the level 
~2x 10~ 5 . The dashed line is drawn for Aa/a = —5.4 x ICT 6 found in the 
MM analysis in 
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this quasar were obtained with the UV- Visual Echelle Spectrograph (UVES) 
installed at the VLT/Kueyen telescope. The VLT/UVES data have a very- 
high signal-to-noise ratio (S/N ~ 90 — 150 per resolution element) and a high 
resolution power of A/AA ~ 55 000 (FWHM ~ 5.5 km s -1 ). The absorption 
systems from this data were analyzed in and j3HJ EH] ■ 

We used (jU)-(|EJ) to estimate the rest-frame equivalent widths W\, the 
line centers A c , and their errors which are listed in Table 1. In this table, 
column (2) presents the oscillator strengths for absorption. The sensitivity 
coefficients Qi and their errors GQ i are calculated from the values qi, a Qi , and 
uj 0ti listed in Table 2 in [2J. 

Our result of the linear regression analysis is shown in Fig. 3 by the 
solid line, while two dotted lines correspond to the 1 a deviations of the 
slope parameter ««. We find z = 1.1490988 ±0.0000016 and Aa/a = (1.1 ± 
1.1) x 10~ 5 (1 a). The uncertainties of the zo and Aa/a values are calculated 
through statistical Monte Carlo simulations assuming that the deviations in 
the Zi and Qi values are equally likely in the intervals [zi — a Zi , Zi + a Zi ] and 
[Qi — dQ^Qi + ctqJ. The result obtained shows that Aa/a = within the 
1 a interval. For comparison, the linear regression for Aa/a = —5.4 x 10 -6 
based on the inhomogeneous sample containing 127 absorption systems j2] is 
shown by the dashed line in Fig. 3. Although the sign of Aa/a is opposite, 
the value of Aa/a = —5.4 x 10~ 6 is consistent with our data within the 2 o 
interval. The accuracy of the regression MM analysis may be increased if 
a few homogeneous samples are combined. The regression analysis is easily 
generalized for this case (see, e.g., [4*3]). 

3 Constraints on the proton-to-electron mass 
ratio 

The proton-to-electron mass ratio /x can be estimated from high redshift 
molecular hydrogen systems. With some modifications, such measurements 
were performed for the z = 2.811 H 2 -bearing cloud toward the quasar PKS 
0528-250 in [H] and |Sj setting the limit of \Afi/fi\ < 1.8 x 10~ 4 (la). Later 
on higher resolution spectra of this QSO obtained with the VLT/UVES 45 
revealed that the H_2 profiles at z = 2.811 are complex and consist of at 
least two subcomponents with Av ~ 10 km s" 1 . Lines arising from the low 
rotational levels with J = and J = 1 show a red-side asymmetry - contrary 
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Figure 4: H 2 absorption features identified at z — 3.025 toward the quasar 
Q 0347-3819 (normalized intensities are shown by dots with la error bars). 
The zero radial velocity is fixed at z H2 = 3.024895. Smooth curves are the 
synthetic H 2 profiles found by the least-squares procedure. The mean S/N 
values per resolution element are indicated 
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to those arising from J = 2 and J = 3 where a blue-side asymmetry is clearly 
seen. This complexity makes impossible higher accuracy measurements of 
A /j,/ ' n in this system. Again, as in the case with the Aa/a estimations, one 
needs a single component narrow H2 absorption-line system to restrict the 
observed wavelengths within the uncertainty interval of a few mA. 

Such system has been found at z — 3.025 toward the quasar Q0347- 
3819 [46 . The absorption spectrum of this QSO was obtained with FWHM ~ 
7 km s _1 and S/N ~ 15 — 30 per resolution element at the VLT/UVES. 15 
unblended H2 lines (selected from ~ 80 identified H2 transitions) which pro- 
vide the most accurate line center measurements are shown in Fig. 4. 

For a given H 2 line, the number of pixels involved in the analysis corre- 
sponds to the number of points in the line profile shown in Fig. 4. To measure 
the line centers we used a method that matches the observed profiles with the 
synthesized ones to estimate a set of model parameters (for details, see |30j). 
The results obtained are presented in Table 2. The measured individual z 
values and their standard deviations are listed in column (5). Using these 
redshifts we can constrain possible changes of \i. 

The rest-frame wavelengths, Ao,i, of different electronic- vibrational-rotational 
transitions of H2 depend in a different way on the reduced mass of the 
molecule. The sensitivity of Ao,i to variation of \i is given by the sensitivity 
coefficient /Q defined as 

fi d\ i 

iCi = - — . (15) 

AQ i a/j, 



Following the procedure described in 04] we calculate the coefficients 
listed in column (6) of Table 2. Their errors are obtained from coefficients 
Y mn (presented in Table 1 in |Sj) using the method of error propagation (Y mn 
values are considered to be accurate to k decimal places and their rounding 
errors are set to 0.5 x 10 _fc ). 

In linear approximation, (JT5j) can be re-written in the form 

= 1 + (/C t - K S ) A/V/i (16) 

or 

Zi = z + {Jd - K) , (17) 
where X ijZ and A Jj2 are the line centers measured in a quasar spectrum, and 

^ = (l + ^o)—, (18) 
A 4 
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Table 2: H2 lines at z = 3.025 toward Q 0347-3819 and sensitivity coefficients 
K 



J Line 


Ag, A 


A c , A 


z 


K 


1... W3-0Q 


947.4218(5) 


3813.2653(41) 


3.024887(5) 


0.0217427(8) 


L1-0P 


1094.0522(51) 


4403.4575(77) 


3.02491(2) 


-0.0023282(1) 


L3-0R 


1063.4603(1) 


4280.3051(30) 


3.024885(3) 


0.0112526(4) 


L3-0P 


1064.6056(5) 


4284.9249(33) 


3.024894(4) 


0.0102682(4) 


L7-0R 


1013.4412(20) 


4078.9785(25) 


3.024898(9) 


0.03050(1) 


L10-0P 


982.8340(6) 


3955.8049(51) 


3.024896(6) 


0.04054(6) 


2... W0-0R 


1009.0233(7) 


4061.2194(61) 


3.024901(7) 


-0.0050567(7) 


WO-OQ 


1010.9380(1) 


4068.9088(53) 


3.024885(6) 


-0.0068461(2) 


W1-0Q 


987.9744(20) 


3976.4943(49) 


3.02490(1) 


0.0039207(2) 


L3-0R 


1064.9935(9) 


4286.4755(98) 


3.02488(1) 


0.0097740(5) 


L4-0R 


1051.4981(4) 


4232.1793(85) 


3.024904(8) 


0.015220(1) 


L5-0R 


1038.6855(32) 


4180.6048(70) 


3.02490(1) 


0.020209(3) 


3... L4-0R 


1053.9770(11) 


4242.1419(28) 


3.024890(5) 


0.012837(2) 


L6-0R 


1028.9832(16) 


4141.5758(31) 


3.024921(7) 


0.022332(7) 


L12-0R 


967.6752(21) 


3894.7974(39) 


3.02490(1) 


0.0440(2) 



"Listed wavelengths are from |4T1 |4*H] . 
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with z and K, being the mean redshift and the mean sensitivity coefficient, 
respectively. 

The linear regression (JT7jl calculated for the complete sample of the H2 
lines from Table 2 provides (A/i//z) 7=1+2+3 = (5.0 ±3.2) x 1CT 5 . However, if 
we consider H2 transitions from individual J levels, then the weighted mean 
redshifts reveal a gradual shift in the radial velocity for features arising from 
progressively higher rotational levels : z(J=l) = 3.024890(2), z(J=2) = 
3.024895(3), and z(J=3) = 3.024904(4). The H 2 lines with changing profiles 
and small velocity shifts as J increases were also observed in our Galaxy in 
the direction of ( Ori A. At present these H2 lines are interpreted as being 
formed in different zones of a postshock gas • 

Although the z(J=l) and z(J=2) values are consistent within 1 a in- 
tervals, the difference between z(J—l) and z(J—3) is essential and equals 
1.0 ± 0.3 km s -1 . This small change in the radial velocity with increasing J 
may mimic a shift in Afi/fi. If we exclude from the regression analysis the 
rotational levels with J = 3, then (A/i//z)j =1+2 = (2.1 ± 3.6) x 10~ 5 . The 
use of lines arising from the same rotational levels would be more reason- 
able to estimate Afi/fi, but in our case the sample size is rather small and 
we have to combine the J = 1 and J = 2 levels to increase statistics. The 
linear regression is shown by the solid line in Fig. 5 while two dashed lines 
correspond to the 1 a deviations of the slope parameter k^. 

Thus, the 1 o confidence interval to A/i//i is —1.5 x 10~ 5 < Afi/^i < 
5.7 x 10~ 5 . For a cosmological model with Qm = 0.3, — 0.7, and H = 72 
km s" 1 Mpc^ 1 , the look-back time for z = 3.025 is 11.2 Gyr jJU]. This leads 
to the restriction \fi/ n\ < 5 x 10~ 15 yr _1 on the variation rate of /i. 

4 Conclusions and future prospects 

We have considered modern experimental investigations of the hypothetical 
variations of the fundamental physical constants from absorption-line spectra 
of distant quasars. The subject was limited to the properties of the two 
dimensionless quantities - the fine-structure constant a and the proton-to- 
electron mass ratio /i. Our main conclusions follows: 

• The error a/Aa/a) — 1-2 x 10~ 6 claimed in |2j is an order of magnitude 
smaller compared to what can be provided by the wavelength measure- 
ments of isolated and unsaturated QSO absorption lines. 
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AyLt/yLi= (2.1±3.6)X10- 5 
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-0.01 0.01 0.02 0.03 0.04 0.05 

Sensitivity Coefficient, K 



Figure 5: Relation between the redshift values z,- L calculated for individual 
H 2 lines shown in Fig. 4 and their sensitivity coefficients /Q. The solid line 
is the linear regression Z{ = zq + ^(/Q — K,). The dotted lines representing 
the 1 a deviations from the optimal slope demonstrate that A///// = at 
the level ~ 5 x 10~ 5 . 
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• A modification of the MM method (regression MM method) which 
makes it more reliable is presented. This modified MM method applied 
to the homogeneous sample of the isolated and unsaturated Fell lines 
yields Aa/a = (1.1 ± 1.1) x 10~ 5 . The uncertainty of this value is in 
agreement with the available accuracy of the wavelength measurements. 

• If A/i/fi = IZAa/a, then the requirement of consistency between 
Aa/a = (1.1 ±1.1) x 1(T 5 at z = 1.149 and Afi/fi = (2.1 ±3.6) x 1(T 5 
at z — 3.025 rules out very large values of \1Z\ > 6, which is in line with 
the results obtained in 123] . 

As stated above, the accuracy of ax ~ 0.2 mA is required to achieve a 
level higher than 10 -5 in the Aa/a measurements. This accuracy is ~ 10 
times higher than that available at modern giant telescopes. In 2003, a new 
High Accuracy Radial velocity Planet Searcher (HARPS) spectrograph 51 
started to operate at the 3.6 m telescope (ESO, La Silla, Chile). Observations 
with HARPS may in principle provide the desirable accuracy ~ 0.2 mA. 
With a resolving power of A/AA = 120 000 and a projection of the fiber over 
4 pixels, one HARPS pixel corresponds to ~ 600 m s _1 . For an isolated 
and unsaturated absorption line, the line center may be measured with an 
accuracy of about 1/20 of the pixel size (~ 30 m s _1 ). This corresponds to 
a\ = 0.4 mA at 4000 A. If isotopic and hyperfine structure systematics as 
well as the others discussed in [2] are known, then we may hope to reach a 
level of ~ 10~ 6 in future astronomical measurements of Aa/a. 
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